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ABSTRACT 

Three transition interpolation schemes for use in h or r-refinement have been analyzed in 
terms of accuracy, implementation ease and extendability. They include blending-function 
interpolation, displacement averaging, and strain matching at discrete points along the 
transition edge lines. The results suggest that the choice of matching depends strongly on 
the element formulations, (viz. displacement or assumed strain, etc.) and mesh refinement 
criteria employed, and to a lesser extent the choice of computer architecture (serial vs. 
parallel) and the equation solution procedures. A recommended paring of some of the 
elements with the choice factors is suggested. 
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1. Introduction 


The h-ve rsion and/or r-version of mesh refinements have been steadily gaining popularity 
in recent years. A primary reason for this popularity has been due to the simplicity 
of data handling for subsequent refinement and solution strategies. The error estimate 
criteria employed for these refinements are by and large deduced from the assumption 
that there is no loss of accuracy in the transitional meshes or the transitional meshes 
are sufficiently remote from the regions of rapid stress or strain gradients (Babuska and 
Dorr, 1981; Felippa, 1976; Shephard and Gallagher, 1979; Zienkiewicz and Zhu, 1987). 
Two techniques most frequently adopted for discretizing the transitional regions are the 
displacement matching and the blending function interpolations (Gordon and Hall, 1973; 
Surana, 1980; Gordon, 1970; Cavendish and Hall, 1984). 

From the viewpoint of element formulations, these techniques go well in hand with the 
displacement-based elements as they involve only displacements at the element boundaries. 
Hence, it appears that transition interpolations for other element formulations such as hy- 
brid and assumed strain elements would require the matching of different field quantities. 
A recent numerical study of mesh refinement strategies that utilize displacement matching 
has shown that the displacement-based 9-node element yields a much better mesh refine- 
ment pattern compared to the assumed natural-coordinate strain (A NS) 9-node element 
(Stanley et ah, 1989). This has motivated us to examine various interpolation schemes 
used in transition elements and classify the desirable conjugacy of element formulations 
vs. transition interpolation schemes. 

2. Formulation of Transition Discretization 


Consider two quadrilateral shell elements that have one common transition boundary as 
shown in Fig. 1. Although one could avoid transition nodes in triangular elements, such 
a strategy can lead to different nodal patterns or severe element distortion, which in turn 
deteriorates the accuracy in the transition regions. This is illustrated in Fig. 2. Conceiv- 
ably, the distortion brought about by an /i-refinement can be regularized to some extent 
if an r- refinement is followed up after each /i-refinement. However, such an r-refinement 
follow-up may not contribute to overall accuracy improvements. 

During the present study of transition discretizations, we will impose one important imple- 
mentation requirement: the construction of transition boundary conditions should involve 
only the nodal degrees of freedom (dofs) along the transition element edges. If one views 
the nodes along the transition element edges as connectors for tying element A to elements 
B and C in Fig. 1, there exist three possibilities. These possibilities include blending 
displacement interpolation, displacement averaging, and strain (or stress) matching. We 
will describe in detail the three procedures. 
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4-Node Element and 9-Node Element and 

Transition Node(X) Transition Nodes (X) 

Transition Node(s) Belong to Element A 


Fig. 1 Transition Nodes for quadrilateral elements due to (h, r)-refinements 



Use of Two Transition Nodes 
Allows Accurate Triangular 
Discretization: 


K=( W 2 


Without Transition Nodes 
Union-Jack Pattern Develops, 
Hence Loss of Accuracy 
in General 


Fig. 2 Transition Nodes vs. Union- Jack Pattern for Triangular /i- Refinements 
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2.1 Blending Displacement Interpolation 

This is by fax the simplest procedure to introduce additional nodes as connectors so that 
the generalized degrees of freedom associated with the additional nodes remain linearly 
independent of the remaining ones. For example, when a 4-node element possesses one 
transition edge, one introduces the following blending interpolations: 


N 

u = 

J=1 

where 

iV, =L 1 (r ? )-Q 1 ( 0 > 
N 2 = LM ■ Q 2 (0 

n 3 = Lj(0 ► 

N a = L 2 (f])-L 2 (0 

If* = Li(n) • Q*it) . 


(i) 


( 2 ) 


in which L's and Q ' s are linear and quadratic Lagrange interpolation functions. A similar 
blending interpolation can be constructed for a 9-node element with one or more transition 
edges. 


T| 



Fig. 3 Blended Interpolation of Transition Nodes 

It should be pointed out that while the interpolation procedure is straightforward, the 
resulting element would require a (2 x 3)-integration rule instead of (2 x 2)-integration in 
order to preserve rank-sufficiency. 
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2.2 Displacement Averaging 


Displacement averaging is by fax the simplest procedure. For example, the displacement 
at node 5 can be obtained by an average of 


u 5 = au 2 + (1 - a)u 2 , 


a = 


ti 

ii + ii 


( 3 ) 


where and £ 2 are the side length from node 1 to node 5 , and node 5 to node 2 , respectively. 
If = t-i, then we have 


$ = u 5 - i(uj + u 2 ) = 0 


( 4 ) 


Hence, 1 x 5 can either be eliminated or the constraint equation (4) can be augmented via the 
method of Lagrange multipliers (A) and u and A can be solved as the system unknowns. 

2.3 Strain Matching 

Instead of displacement matching or interpolation, one may impose strain matching as 
a connector requirement. For this case, in conjunction with the present implementation 
requirement, viz., the use of dofs along the transition edges only, we axe to match the 
membrane, transverse shear and bending strains: 




7e 




1 

du^ 


di 

1 

dw^ 

* 

d( 

1 

dd 

A* 

ft, 


( 5 ) 


where is the first fundamental coefficient of the neutral surface of the shell along the 
boundary £, are the tangential, normal and rotational deformation along 

For a 4-node element having one transition edge as shown in Fig. 3 , it turns out the strain 
matching is the same as the displacement averaging as long as the transition node 5 bisects 
nodes 1 and 2 . To examine the strain matching for a 9 -node case, we adopt a flat plate 
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for illustrative purposes as shown in Fig. 4. 


"H 



Fig. 4 Strain Matching Along A 9-Node Edge 

Since there are two transition nodes with (u, w, 0,j)-dofs affected, a totcil of six constraints 
axe required in order to eliminate them. Hence, we must enforce the strain matching of the 
three strain states given by (5) at two locations along the transition edge. Among several 
candidates, it has been determined that the most desirable strain-matching locations are 
the two Barlow points (£ = ±l/\/3). Matching the strains at the two Barlow points yields 
the following constraints: 

f 8 u 4 — (3uj + 6 u 2 — u 3 ) = 0 
8 u 5 — (—ui + 6 u 2 + 3 u 3 ) = 0 
8t>4 — (3tq + 6 u 2 — v 3 ) = 0 
8v 5 — (—vi + 6 u 2 + 3u 3 ) = 0 
804 - (30i + 60 2 - 0 3 ) = 0 

4> = < 80s — (—0i 4- 60 2 + 30 3 ) = 0 ( 6) 

80 4 - (301 4 602 -0 3 ) = O 

80s -(-0i +602+ 30s) = 0 

e 

8w 4 — (3u)i 4- 6u) 2 — w 3 ) + -(03 — 20 2 + 0i ) = 0 

i 

8u> 5 — (— u>i 4 6u>2 + 3w 3 ) — — (0 3 — 20 2 + 0i ) = 0 
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Rewriting (6), one obtains 


Bsd = 0, d= [ui «i t»i 8\ ft, u 5 v 5 w 5 6 5 \ T (7) 

If the displacement averaging is adopted, one replaces the two strain constraint equations 
on the transverse shear strains by 


8w 4 — (3wj + 6 u >2 — w 3) = 0 
8 w 5 — (—w i + 6 u >2 + 3 ws ) = 0 


( 8 ) 


3. Implementation 

There axe two distinct needs for implementing the blending function approach vs the 
method of Lagrange multipliers. We shall describe the two procedures separately. 

3.1 Blended Assumed Natural-Coordinate Strains 

To illustrate the blended assumed natural-coordinate strain (ANS) interpolation, let us 
consider the following 5-noded membrane square element. 


n, y, v 



Figure 5 5-Node Blended- Function Square Plate Element 
From the above figure one readily obtains the following strain expressions: 
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( 9 ) 


To see whether (9a) yields a rank-sufficient stiffness matrix, it is adequate to check 


K f 


= / B^ X EB XX dV 

Jv 


( 10 ) 


where 


Cxx = B zx U, U=L«1, «2, u 3i «4, «5j T (11) 

B xx = [ ^(-1 + j(\ + 2£)Ni{r)), — jN 2 (t]), jN 2 (t)), ( 12 ) 

in which 


WOO =§(!-*?), N 2 ( V ) = ^1 + T,) 


and E is Young’s modulus. Evaluation of (10) for a constant thickness plate ( h ) gives 
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( 13 ) 


which has the required correct rank since there is only one physical rigid mode. Hence, 
(2 x 2)-integration of the transition element involving node 5 will have a correct rank. 
For higher order elements, e.g., 9-node element, similar ANS transition elements can be 
constructed. 


3.2 Augmentation of Displacement Averaging or Strain Matching Condition 

For linear cases, the transition conditions given by (6) or (7) can be augmented as follows: 


<m = <5u r Ku + <$A t Bu + £u t B t A - <5u T f = 0 


so that one obtains 


K 

B 


B Tl 
0 


]{:}-« 


(14) 

(15) 
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Partitioning u further into the transition node u* and the original node Uo, (15) can be 
rearranged as 


K 00 

K ot 

an 


u ° } 


[ f o] 

K to 

K tt 

Bf 

< 

u ‘ 

► = < 

f< 

B„ 

Bj 

0 


i A J 


0 J 


(16) 


So that (15) can be solved directly or by 


hojjUo — fo 


(17) 


where 


K„„ = K, 




K ot Br 1 Bo-Bj’B t - T K to 

T t>—T t 


+ B^B ( - T K u Br 1 B 0 


-1 


(18) 


4. Strain Matching for 9-ANS Transition Element 

As alluded to in Introduction, the ANS family of shell elements are constructed by inter- 
polating the strains along the element edges and certain interior natural-coordinate lines. 
Hence, displacement matching or blended function approach can destroy much of the good 
ANS element attributes. A preferred construction of transition strategies for ANS elements 
is to match the strains along the transition edges at some discrete points. To this end, let 
us consider Fig. 6 in which the element edge A is to match with the element edges B and 
C as shown below. 

T\ 
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The position vector of the particle point P (Fig. 7) on the deformed shell is given by 

r = r° + Cb 3 (19) 


where 

r° = xei + ye 2 + ze 3 , (x y z) = {(X + u) (F + u) (Z + w)} (20) 

in which x, yand z are the deformed neutral shell surface position coordinates, (u, vand w ) 
are the displacements measured in the inertial e-system, £ is the distance of the material 
point P from the shell neutral surface measured in the b-system attached on the deformed 
cross-section of the shell, and the vector b is related to the vector e by 

b = Re (21) 


The linearized (or incremental) displacement vector, U, is obtained from: 

u = u r e + £u (22) 

where the pseudo-rotation vector, u, is related to the shell-surface, pseudo-rotation quan- 

tities, (3 according to 

u=(-tr,(2) £(!)>{*}, P = Tl& T„ (23) 

in which tf 9 (i) represents the i-th row of the transformation matrix, T as , defined by 

s = T sg e (24) 

that are attached to the deformed shell surface; and, Tfc s relates the shell-surface basis 
vectors s, to the inertial basis vectors according to 

b = T fca s (25) 


The covariant natural- coordinate unit vectors are then obtained by: 

1 dr° 1 


a f = 


At d£ At 


- 7 - (x,( ej + y,t e 2 + z € e 3 ) = t^ e 


1 dr° 1 , , xT 

a„ = — -x— = — (x,, ei + y, v e 2 + z,e 3 ) = t, e 


" A v dy A v 


&£ X d-n / V T 

a< = t 7 = (x,c ej + y ?< e 2 + z, c e 3 ) = t c e 

\Ci£ X Gjj I 


(26) 


(27) 

(28) 
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Figure 7 Geometry of Shell Surface 


where the two fundamental shell surface quantities, and A,, are given by: 


2 _ dr^ dr^ 2 dr° dr° 

^ d£ d£ ' v dr) dr) 

For subsequent applications, we express the above relation in a compact form: 


(29) 


a — T n gG — 



Finally, the covariant partial derivatives are given by: 


(30) 


d_ld d_ld did 
Wf. ~ A~tdi dS ~ ~ A,d^ dS( ~ A c d( 


A< = h(Z, Tj)/2 (31) 


where /i(£, T)) is the shell thickness. 

First, the membrane strains, e^, at the two Barlow points along the edge A are given ( 
see Park et al.(1989) for details in the natural- coordinate strain formulation), respectively, 
by: 


c «(&) = ^- t f,-[|( u 3-«i)+6(«3-2u 2 + u 1 ) (32) 

Zb 
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where £& and ^ are the two right and left Barlow points along the edge A. 

On the other hand, the assumed strains at the same location along the edge B and C can 
be derived, respectively, as 


S6 Sfc 


(34) 


&'(&) — IC 1 + 77) e ^'(^) + 1(1 - 7r) e ?«'(^6 


(S) 


(35) 


in which ^f>^(Cb) an d e ^'(^) 8X6 strains at the two Barlow points along the edge C, 
and e®£»(£j) and (£j) axe the strains at the two Barlow points along the edge B. 

Equating e^(^t) and one obtains the following constraint at £(, = l/\/3 on the 

edges A and C: 


cl ■ u 5 + ■ u 3 + cl ■ u 2 + cf • Ui = 0 


(36) 


where c* are given by 




( Cl = (0.5 - 


where the superscript C denotes the quantities pertaining to the edge C. 

Similarly, the membrane strain matching constraint at = — 1/\/3 along the edges A and 
B is obtained by: 



cl ■ u 4 + ■ u 3 + cl - U 2 + cl ■ u, = 0 


(38) 


where Cj are given by 



{ c 4 = (1 - 1/V3)^tg + (1 - y/3)-fctf. 


' 6 


c 3 =(0.5-l/v/3)^t^ 
e2 = + 1 /^)^r t § + j(-5 + 3 \/3)^-t| 

Cl = -(0.5 + l/%/3)^-t^ + i(-3 + 5/>/5)^rt €; + ±(1 + v/3)^-t| 


( 39 ) 


It should be noted that the two bending strain constraints at the two Barlow points 
axe obtained by replacing u by u in the preceding membrane matching conditions. The 
remaining two strain matching conditions, viz., the transverse shear strains, 


1 <9u 


T'€ = t C^-^7+ t [» 


C A, 


lead to the following two constraints: 


(40) 


dg • U 5 + d, • u 3 + d^ • U 2 + dj • Uj + e^ • U 5 + ef • U 3 + e^ • U 2 + ef • Ui — 0 (41) 


,T 


-,r 


u 4 + d^ • u 3 + d^ • u 2 + d l • U! + e* • u 4 + e J ■ u 3 + ■ u 2 + ef • u x = 0 (42) 


-T ~ 


-T 


where dj and dj are obtained by replacing with t^, etc., into c, and Cj, respectively, 
and ej and ej are given by 


e 5 = (1 - 1/V3)tg + ±(>/5 - l)tg 
e 3 = - 5(1 + n/ 3)% + + £(\/3 ~ 2 ) t ^ 

e 2 = 1% + |(3 - 2\/3)t^ + |t^ 
l ei = -|(1 - 


' e 4 = i(v/3 - l)tf ; + (1 - l/>/5)tf ; 

«3 = ~J(1 ~ V^)t^ 

e 2 = + |(3 - 2\/3)t^ 

, ci = -|(1 + V3)t 6 + |(\/3 - 2)t^ + 


(43) 


(44) 


Finally, the displacement u n and u n , namely, those perpendicular to and t^ ( t„ 
x t^), can be obtained by 


gif • u 5 + gf • u 3 + gif • u 2 + gf • Ui =0 


(45) 
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where 


gif • “4 + g;T • U 3 + • U 2 + gf • U, = 0 


(46) 



The interpolation functions for u n at the 


- N 

g4 


' 8 ' 

g 3 

g2 

> = t„(f 6 ) < 

1 

-6 

k gl J 




Barlow points axe the same as for u n . 


(47) 


Combining the six strain matching constraints and the four displacement interpolation 
conditions together, one obtains the following equation: 


*s = B ans { q } 

where q t = (u t -, u t ) T and B ans is given by 


B ans 


K 

0 

T 

c 2 

0 

C T 

c 3 

0 

0 

0 

c T 

c 5 

0 

0 

T 

C 1 

0 

T 

C 2 

0 

T 

c 3 

0 

0 

0 

T 

C 5 

dl 

T 

e l 

a t 

d 2 

e T 

e 2 

a t 

d 3 

T 

e 3 

0 

0 


c T 

c 5 

si 

0 

T 

g2 

0 

T 

g 3 

0 

0 

0 

g 5 T 

0 

0 

T 

gl 

0 

T 

g2 

0 

T 

S 3 

0 

0 

0 

T 

si 

-T 

C 1 

0 

C T 

c 2 

0 

c T 

c 3 

0 

-T 

C 4 

0 

0 

0 

0 

-T 

C 1 

0 

-T 

c 2 

0 

-T 

c i 

0 

C T 

l 4 

0 

0 

df 

-T 

e l 

a t 

d 2 

-T 

e 2 

a t 

^3 

-T 

e 3 

a t 

d 4 

-T 

e 5 

0 

0 

-T 

gi 

0 

-T 

g2 

0 

-T 

si 

0 

-T 

si 

0 

0 

0 

. 0 

-T 

gl 

0 

-T 

g2 

0 

-T 

si 

0 

-T 

si 

0 

0 


(48) 


(49) 


It can be shown that the above expression reduces to the fiat, straight edge case given by 
(6). It should also be noted that u introduced in (23) must be replaced by (fa, fa) in (48) 
in actual implementation in order to avoid singularity. 


5. Discussions 

The present study has examined three possible ways for discretizing the mesh transi- 
tion regions. It has been shown that for four-noded elements the displacement averaging 
and the strain matching yield the same transition discretization expressions. For plates 
with straight edges discretized by nine-noded elements, the Barlow-point displacement and 
strain matching becomes the same except for the transverse shear strains. 

The transition interpolation schemes for curved shell boundaries for the nine-noded ANS- 
shell element (Park and Stanley, 1895; Stanley, 1985; Park et al., 1989) has been carried 
in detail for implementation in h-version mesh refinement. It is noted that the transition- 
boundary strain matching directly affects the resulting assumed natural-coordinate strains 
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employed in the formulation of ANS-shell element generation. Such explicit enforcement 
of strain-matching conditions along the element edges is considered a distinct advantage 
of the ANS-shell element. 
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